Energy Barriers to Motion of Flux Lines in Random Media 
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We propose algorithms for determining both lower and upper bounds for the energy barriers 
encountered by a flux line in moving through a two-dimensional random potential. Analytical 
arguments, supported by numerical simulations, suggest that these bounds scale with the length t 
of the line as t 1 ^ 3 and t 1 ^ 3 \f\at, respectively. This provides the first confirmation of the hypothesis 
that barriers have the same scaling as the fluctuation in the free energy. 



Glassy states are characterized by a complex energy 
landscape with many (metastable) free energy minima 
[0,1). Some commonly encountered examples are spin 
glasses, pinned flux lines in superconductors, polymers, 
domain walls in random field and random bond Ising 
models, and many more. Typically, the fluctuations be- 
tween free energy minima in these systems (either in dif- 
ferent realizations, or in the same random system) scale 
with observation size L as L e . On the other hand, dy- 
namic response of the system is limited by the barriers in 
free energy encountered in crossing from one minimum to 
another. It is likely that the scale of these barriers also 
grows with size as . The simplest assumption is that 
the two energy scales are comparable, i.e. tp = 9. How- 
ever, it is also quite possible that the heights of the ridges 
in the random energy landscape scale differently from 
those of the valleys that they separate, with tp > 9. Yet 
another scenario is that transport occurs mainly along 
a percolating channel of exceptionally low energy valleys 
with ip < 9. In this paper we examine a specific glassy 
system, a flux line moving in a two dimensional random 
potential, for which we demonstrate ip = 9. 

The system we study is inspired by measurements 
of nonlinear current-voltage (J, V) characteristics, V oc 
exp(— const of disordered superconductors in a 

magnetic field || . There is emerging consensus that such 
behavior is best described in terms of low-temperature 
glass phases MM: In the weak current regime, the dy- 
namical properties are dominated by activated pro- 
cesses corresponding to bundles of flux lines overcoming 
pinning barriers. The principal difference from the classi- 
cal picture of flux line motion is the power-law growth 
of barrier energies Eb oc L^, with bundle size L. The 
latter in turn diverges as the superconducting current 
driving flux flow decreases, leading to nonlinear (/, V) 
dependence quoted above. As accurate values of critical 
exponents characterizing vortex glass phases (including 
ip) are presently not available, direct experimental con- 
firmation or refutation of the various theoretical models 
is difficult. The point of view taken in many of the pio- 
neering papers |8|,|| on the subject is that ip = 9. Here we 



shall attempt to put this assumption on a firmer basis. 

We shall examine the configurations of a single flux line 
(FL) in a random potential landscape. Equilibrium prop- 
erties of this system have been extensively studied in the 
context of directed polymers in random media (DPRM) 
. It is known that the FL is pinned by impurities into 
a glassy state. Furthermore, by using a transfer matrix 
method, properties of this state can be probed numer- 
ically in polynomial time in the line length t. In two 
dimensions, e.g. for a FL trapped between two copper 
oxide planes of a high-T c superconductor, some analyti- 
cal information is also available. For example, the fluc- 
tuations in the free energy at finite temperature scale 
as i 1 / 3 . Since the scaling behavior of the pinned FL is 
governed by a zero-temperature fixed point g), energy 
fluctuations scale in the same way. The availability of 
such results and techniques make this system ideal for 
investigation of barrier energies. 

The precise description of the model is as follows: The 
FL is discretized to lie on the bonds of a square lattice, 
directed along its diagonal. Each segment of the line can 
proceed along one of two directions, leading to a total of 
2* configurations after t steps. These configurations are 
labelled by the set of integers {x(t)} for r = 0, 1 , • • • , t, 
giving the transverse coordinate of the line at each step 
(clearly constrained such that x(t + 1) = x{r) ± 1). To 
each bond on the lattice is assigned a (quenched) ran- 
dom energy equally distributed between and 1. The 
energy of each configuration is the sum of all random 
bond energies on the line. Without loss of general- 
ity, we set x(0) = 0. For each endpoint (t,x) with 
x = —t, — t + 2, • • • , t, there is a configuration of mini- 
mal energy E m i n (x\t) which can be obtained numerically 
in a time of order t 2 . It is known that for \x\ < x c oc t 2 / 3 
the function E m i n (x\t) behaves as a random walk and is 
thus asymptotically Gaussian distributed |^[|. We next 
examine the energy barrier that has to be overcome when 
the line is moved from an initial minimal energy configu- 
ration between (0,0) and (t, —Xf) to a final one between 
(0,0) and (t,+x f ). 

The only elementary move allowed is flipping a kink 
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along the line from one side to the other (except at 
the end point). Thus the point (t,x) can be shifted to 
(t, x ± 2). Each route from the initial to the final con- 
figuration is obtained by a sequence of such elementary 
moves. For each sequence, there is an intermediate con- 
figuration of maximum energy, and a barrier which is the 
difference between this maximum and the initial energy. 
In a system at temperature T, the probability that the 
FL chooses a sequence which crosses a barrier of height 
Eb is proportional to exp(— Eb/T), multiplied by the 
number of such sequences. We assume that, as is the 
case for the equilibrium DPRM, the "entropic" factor of 
the number of paths does not modify scaling behavior. 
Thus at sufficiently low temperatures the FL chooses the 
optimal sequence which has to overcome the least energy, 
and the overall barrier is the minimum of barrier energies 
of all sequences. 

Since the number of elementary moves scales roughly 
as the area between the initial and final lines, the num- 
ber of possible sequences grows as t xt . This exponential 
growth makes it practically impossible to find the bar- 
rier by examining all possible sequences, hampering a 
systematic examination of barrier energies. Rather than 
finding the true barrier energy, we proceed by placing 
upper and lower bounds on it. The lower bound was 
given in ref. |10|, and scales as t 1 / 3 . In this paper, we 



present an algorithm for obtaining an upper bound. An- 
alytical arguments, suggest that this upper bound grows 
as t 1 / 3 v / hi7, thus establishing ijj = 1/3. Since these argu- 
ments do not constitute a rigorous proof, we verify their 
validity by numerical simulations. Computer time and 
memory requirements for the construction of this upper 
bound are happily polynomial in t. 

A lower bound for the barrier energy is obtained as 
follows PQ |: Since the endpoint of the path has to 
visit all sites (t,x) with |x| < Xf, and since the en- 
ergy of any path ending at (t, x) is at least as large 
as E m i n (x\t), the barrier energy cannot be smaller than 
m&x[E mm (x\t) - E min {-xf\t)] for x <E [— Xf,x/]. When 
Xf is sufficiently small, the probability distribution of this 
lower bound is identical to that of the maximal deviation 
of a random walk of length Xf. The latter is a Gaussian 
distribution with a mean value oc ^Tx] and a variance 

<x Xf. This growth saturates for xf of the order of i 2 / 3 , 
leading to the scaling behaviors, 



To obtain an upper bound for the barrier, we specify 
an explicit algorithm for moving the line from its initial 
to its final configuration. This is achieved by finding a 
sequence of intermediate steps. It is certainly advanta- 
geous to keep the intermediate paths as close to minimal 
configurations as possible. We first attempt to move the 
path in steps from a minimal configuration with endpoint 
at (£, x) to one with endpoint at (t, x + 2), starting at 
(t, —Xf) and ending at (t 7 Xf). At each step, we obtain a 
local barrier path which separates two neighboring min- 
imal paths. The overall barrier is of course the one with 
the highest energy. While it may occasionally be pos- 
sible to go from one optimal path to a neighboring one 
in a single elementary move (as defined above), this is 
generally not the case. Minimal paths with neighboring 
endpoints may be quite far apart at coordinates r < t. 
The reason is simple: suppose the random potential has 
a large positive fluctuation, a "mountain" . The minimal 
energy paths will then circumvent this region by going to 
its right or left. The last path going to the left and the 
first one going to the right have almost the same energy 
(these energies are strictly equal in the continuum limit). 
They form a loop which can be quite large and is likely 
to enclose the barrier when both paths separate already 
at small r. Such loops have been conjectured to 
play an important role in the low-temperature dynamics 
of the DPRM. Since the transverse fluctuations of a min- 
imal path of length t grow as t 2 / 3 , we expect the lateral 
size of these loops to also be of this order. 
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FIG. 1. Scaling functions /i(y), f 2 (y), gi{y), and g 2 (y) (see 
eqs.(l), (4), and (5)) for the mean and variance of the lower 
and upper bounds; averaged over 2000 realizations of random- 
ness, for t = 256 (solid), t = 512 (dotted), t = 1024 (dashed), 
t = 2048 (long dashed), and t = 4096 (dot-dashed). 



for the lower bound and its variance. The functions fi(y) 
and /2(y) are proportional to ^Jy and y for small y, re- 
spectively, and go to a constant for y = 0(1). Our simu- 
lation results for systems with t = 256, 512, 1024, 2048, 
and 4096 confirm this expectation. Fig. |l| shows the scal- 
ing functions fi{y) and f 2 {y) for different t, and the col- 
lapse is quite satisfactory. However, the initial growth 
oc y/xj, is not clearly seen at these sizes. 



The algorithm for moving a line of length t = 2" from 
a minimal configuration {xi(t)} to another one { 3:2(1")}, 
with x 2 (t) — xi(t) + 2 is as follows: If 2:2(1") < xi(t) + 2 
for all r, we can choose a sequence of elementary moves 
such that at most two bonds of the line are not on one 
or the other minimal path, leading to a barrier of order 
1 between the two. If X2(t) > x\{t) + 2 for some r, the 
two paths enclose a loop. We then consider the midway 
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points (t/2,x) with x%(t/2) < x < X2(t/2). For each of 
these points, we find two minimal segments of length t/2 
connecting on one side to (0, 0) and on the other to either 
(t,xi(t)) or (t,X2(t)) The two segments form an al- 
most minimal path of length t, constrained to go through 
the point (t/2,x). We next move the line {xi(t)} step- 
wise through this sequence of almost minimal paths. At 
each step we first attempt to move the upper segment 
and then the lower one . The prescription for moving 
these segments of length t/2 is exactly the same as for 
paths of length t: If the distance between two consecutive 
configurations is larger than 2 for some r £ [0, t/2], we 
consider the points at (t/A, x) in between the two, and 
find minimal paths of length t/A connecting them to the 
initial and final points. Next we attempt to move seg- 
ments of length t/2 by repeatedly moving line portions 
of length t/A. In some cases, when the energy barrier 
is high, it is necessary to proceed with this construction 
until the cutoff scale t/2 n — 1 is reached. Thus, at each 
intermediate configuration, the line is composed of one 
minimal segment of length t/2, one of length t/A, etc; 
ending with two smallest pieces of length t/2 m (equal to 
1 in the worst case). The barrier paths resulting from 
this construction, and the minimal paths separated by 
them, are shown in Fig. ^. 




FIG. 2. Minimal paths of length t — 256 to endpoints 
between x — —96 and x — +96 (solid), and the barrier paths 
between them (dotted). 

We now estimate the barrier energy resulting from the 
above construction. Each intermediate path is composed 
of segments of minimal paths with constrained endpoints, 
and we would like to find the probability distribution 
for the highest energy. Constraining the endpoints of a 
minimal path of length r typically increases its energy 
by E_{t) cx t 1 / 3 . A subset of these intermediate paths 
(those that cross the largest mountains) have constraints 
imposed on segments of length t, t/2, t/A, and all the 



way down to unity. The number of paths in this subset 
(henceforth referred to as candidate barriers) grows as 
N c (t) cx t a , with l<a<l + 2/3. The lower limit comes 
from noting that at each bisection of a loop, several new 
large loops are generated, at least one in the upper and 
one in the lower half of the parent loop, thus N c > t. The 
upper limit comes from the total number of intermediate 
steps that grows as txf. The energy of each candidate 
barrier path is obtained in a manner similar to that of 
the lower bound: Instead of finding the maximum of a 
random walk of length Xf cx i 2 / 3 , we now have to exam- 
ine the sum of the maxima for a sequence of shorter and 
shorter random walks added together. The mean value 
of this sum is related to the convergent series, 

(E c (t)) = {E-(t)+E- (t/2) +£L(t/4) + .-.) = 
(£L(t)) (l + 2- 1 / 3 + 2~ 2 / 3 + •••)+ A ~ 

(£-(*)) (1-2- 1 / 3 ) 1 +A = AM...(E_(t)) +A. (2) 

The constant A in eq.(|^) accounts for the breakdown of 
the scaling form of the energy increase for small loops. 
The mean angle of the smallest loops (of size t m = t/2 m ) 
approaches the 45° limit; their mean energy growing as 
0.5t m . For the larger loops, the angle tm 3 /t m is small and 
the energy is 0.23t m . A finite value of m acts as a cutoff 
separating the two limits. The energy difference per unit 
length between small and large paths then leads to the 
additive constant A (of the order of unity) in eq. (|J) . 

The barrier energy is the maximum of the N c (t) ener- 
gies of all candidate barriers. To find its characteristics, 
we need the whole probability distribution for the en- 
ergy E c (t). Since E c is the sum of energies coming from 
its minimal segments, the simplest assumption is to re- 
gard the segment energies as independent, approximately 
gaussian, random variables. We then conclude that E c (t) 
is also gaussian distributed with a variance, 

var (E c (t)) = var (-EL(t)) + var (£L(t/2)) H 

~ 2.70...var(£L(t)) cx t 2/3 . (3) 

Since the different segments are in fact constructed 
through a specific recursive procedure, their indepen- 
dence cannot be justified. Thus the statements of the 
gaussian nature of E c (t), and the variance in eq.(|^), 
should be regarded as plausible assumptions that appear 
to be supported by the numerical simulations. 

It can be checked easily that (for large N), the maxi- 
mum of N independent gaussian variables of mean a and 
variance a 2 , is a gaussian of mean a + a\ // 2\nN and vari- 
ance <r 2 /(2 In TV). Since the N c (t) candidate barriers have 
large segments in common, their energies are not inde- 
pendent. We can approximately take this into account 
by assuming a subset of them as independent, leading to 
N cx t a for some a' < a. We thus obtain the following 
estimates for the mean upper bound in barrier energy, 

(E+{x,t)) = (E c (x,t)) + ^2\nNv&rE c {x,t) 
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and its variance, 
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The functions <?i(y) and 32(2/) are proportional to ^Jy and 
y, respectively, for small y, constant for large y, and in 
general different from that of the lower bound. 

Our numerical simulations indeed confirm the above 
scaling forms. The scaling functions gi(y) and 52(2/) are 
plotted in Fig. [I] for different values of t, after av erag- 
ing over 2000 realizations of randomness. The ^/ln(t) 
factors are essential, as a comparable collapse is not ob- 
tained without them. In fact the best fit to < E + (t) > 
is obtained by including the correction to scaling factor 
4.85 < E_(t) >, and with /3 = 1. The numerics therefore 
support the neglect of correlations, and the assumption 
of a gaussian distributed E c (t). As in the lower bound, 
the initial scaling cx y/xj is not clearly seen for the sizes 
studied. Since the leading power for the scaling of the 
lower and upper bounds is identical, we conclude that 
the barrier energies also grow as t 1 ^ 3 . (It remains to be 
seen if the logarithmic corrections are truly present, or 
merely an artifact of our algorithm.) 

We now return to the original question of the response 
of a flux line to an external force, which in the context of 
superconductivity is proportional to the supercurrent I. 
The standard argument (6|4j assumes an exponential de- 
pendence of the net velocity on the typical barrier height. 
However, it is quite possible that the overall response of 
the system is determined by the largest, rather than the 
typical barriers. If so, knowledge of the probability dis- 
tribution of energy barriers is important. For example, it 
may be more appropriate to average the waiting time for 
the activated jumps over barriers, r cx exp(i?B/T). As- 
suming a gaussian decay in the tail of the barrier heights, 
the latter average, [r] a „ cx J dEBe EB / T P(Eb), is domi- 
nated by energies E* B oc with ijj = 2ip = 2/3. More 
detailed considerations of such issues will be taken up in 
future publications jl3[ |. 

In conclusion, for the simple example of a DPRM, we 
have shown that fluctuations in the minima of the energy 
landscape, and the barriers between them, both scale 
with the length of the line as i 1 / 3 . It remains to be 
seen if the upper bound can be further improved upon, 
and placed on a more firm analytical basis. These results 
provide a glimpse into the complexity of the free energy 
landscape of glassy systems. 
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Note: After submission of the manuscript we learned 
of similar methods being developed by A. Middleton. 
One of us (BD) succeeded in generalizing our algorithm 
to a FL in 3 dimensions, again confirming ip — 9 (see 

0). 
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